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i 1 ABSTRACT 

rf\ Context. Our understanding of the formation and evolution of planetary systems is still fragmentary becuase most of the current data 
provide limited information about the orbital structure and dynamics of these systems. The knowledge of the orbital properties for 

i-G a variety of systems and at different ages yields information on planet migration and on star-planet tidal interaction mechanisms. 
C*H Aims. In this context, a long-term, multi-purpose, observational programme has started with HARPS-N at TNG and aims to 
JL characterise the global architectural properties of exoplanetary systems. The goal of this first paper is to fully characterise the 
orbital properties of the transiting system Qatar-1 as well as the physical properties of the star and the planet. 
Methods. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter- 
McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric- 
transit light-curves were analysed and a spectroscopic characterisation of the host star atmospheric parameters was performed based 
on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). 

£SJ Results. We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment 
^ of the system; this information, combined with the spectroscopic determination of the host star properties (rotation, T e g, log g, 

l/~} metallicity), allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution 

f*"*) for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of A = —8.4 ± 7.1 deg. 

C"^*) The planet, with a mass of 1.33 ± 0.05 M.j, is found to be significantly more massive than previously reported. The host star is 

(""*) confirmed to be metal-rich ([Fe/H]= 0.20±0.10) and slowly rotating (vsin/ = 1.7 ± 0.3kms _1 ), though nroderately active, as 
indicated by the strong chronrospheric emission in the Call H&K line cores (logR^K ~ —4.60). 

Conclusions. We find that the system is well aligned and fits well within the general A versus T e s trend. We can definitely rule 
out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present 
orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions. 
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1. Introduction 



The study of extrasolar planets and the determination 

Send offprint requests to: E. Covino of their observational properties have made remarkable 
Based on observations collected at the Italian Telescopio gress over the ast decade . The surve s conducted SO 

Nazionale Galileo ( 1NG), operated on the island of La Palma by i .,, ,, r , , , . ,. ,. , 

4-u -c j •- n Vi A i-i • r u. tat aw it iv at • l f ar with the most succcsslul techniques (i.e. radial veoci- 

the lundacion Galileo Galilei oi the 1JNA1 (Istituto Nazionale . . . . si 111 

di Astrofisica) at the Spanish Observatorio del Roque de los ties and planetary transits photometry) have revealed plan- 

Muchachos of the Institute de Astrofisica de Canarias, in ets spanning an unexpectedly wide range of orbital prop- 

the frame of the programme Global Architecture of Planetary erties. Particularly surprising is the detection of close-in 

Systems (GAPS). 

** Also based on observations collected at Asiago Observatory, strasbg.fr (1 30.79.128.5) or via |http://cdsweb.u-strasbg.fr/cg>| 

and Calar Alto Observatory. Table[3] is only available in elec- bin/qcat?J/A+A/| 

tronic form at the CDS via anonymous ftp to cdsarc.u- Correspondence to: covino@oacn.inaf.it 



Covino et al.: GAPS. I. RM effect in Qatar- 1 



giant planets with orbital periods as short as one day, the 
so-called hot Jupiters. 

Today ground-based and space-borne photometric sur- 
veys are yielding cr ucial information on transiting planets 
( Batalha et al. 2013 ). However, the radial velocity technique 
is still not only of essential value (e.g., to confirm transit- 
ing planet candidates and determine their actual masses), 
but is now beginning to extend searches to still unexplored 
ranges of planet and host star properties, thanks to the 
improved sensitivity and stability of the new generation 
of high-resolution spectrographs, pioneered by the HARPS 
instrument on the 3.6m ESO telescope, and now followed 
by the newly built high-resolution spectrograph HARPS-N, 
recently come into operation on the Telescopio Nazionale 
Galileo, TNG ( |Cosentino etllpoTS ). 

One of the relevant open questions in the exoplanetary 
field concerns the characterisation of the architectural prop- 
erties of extrasolar planets and their possible dependence 
on the physical properties of the parent stars. For example, 
known systems with giant planets are usually not followed- 
up to search for additional low-mass companions that might 
exist at a range of separations: as a consequence, we still 
lack information about the frequency of solar-system-likc 
systems. Properties such as frequency and orbital charac- 
teristics of exoplanets are expe cted to depend upon stellar 



et al. 



properties, such as metallicity (ISozzetti et aLl|2009 Santos 
2011 Mortier et a l. 2013J and mass (Johnson et al 



2010 



Bonfils~et al.||2013D, as well as on the host star's en 



vironment QDesidera fc Barbieri|2007 Pasquini et al.|20"l2 



Quinn et al.||2012| |Roell et al.||20T2p Furthermore, the rel- 



ative role of different mechanisms of the time evolution of 
planetary system architecture is still rather uncertain. 

Crucial information on the mechanisms governing the 
evolution of planetary systems can be obtained by deter- 
mining the relative orientation of the host-star spin axis, 
which is usually regarded as a relic of the angular mo- 
mentum of the protostellar accretion disc, and the orbital 
axes of the planets. Of the two main mechanisms invoked 
to explain inward migration of giant planets from their 
original location, disc-planet interactions tend to preserve 
the initial spin-orbit alignment as well as orbit circular- 
ity, wherea s dynamical interactions, such as planet-planet 



scattering (Papaloizou & Terquem 2006 Chatterjee et al. 



2008) or Kozai resonanc e due to the presence of an off 



plane massive perturber ( Fabrycky fc Tremaine| [2007), are 
expected to alter both the inclination and eccentricity of 
the orbit. Planet-planet scattering must occur after the 
disc dissipated, otherwise the disc in teraction would drag 
the planet back on the m edian plane ([Bitsch fc KleypOll 



Marzari fc Nelson 2009[). A lthough, as suggested recently 



by |Teyssandier et al. ( |2013 ), for a sufficiently high initial 
inclination of the orbital plane, even a disc could lead to 
Kozai cycles, and hence to spin-orbit misalignment. The 
orbital eccentricity can become very high during Kozai cy- 
cles with consequent decrease of periastron distance and 
set-in of star-planet tidal interactions unt il the orbit even- 
tually becomes very tig ht and circularised ( Katz et al.|2011 
Lithwick k Naoz|2011 ) 



Determining the orientation of a planet orbital axis with 
respect to the stellar spin axis is therefore a way to as- 
sess the relative importance of the two migration mecha- 
nisms. This can be accomplished through the observation 
of the Rossiter-McLaughlin (R M) effect, well known from 
the study of eclipsing binaries (Rossiter]|1924 McLaughlin 



1924). The RM effect is an anomaly in the radial ve- 
locity (RV) curve occurring during a planetary transit, 
whose shape yields information on the sky-projected an- 
gle A between the star spin axis and the planet orbital axis. 
Therefore, measuring the RM effect provides a unique ob- 
servational constra int to the actual spin-orbit misalignment 
( Winn et al.|2010" ). In some cases, the analysis of starspots 
can also provide a good determination of the sky-projected 
obliquity A and, in the most favoura ble ones, can even yield 



the true obliquity (V>) of the orbit (|Nutzman et al. 
Tregloan-Reed et al.|[20T3 | 



2011 



So far, A has been measured for a growing number of 
transiting planet JH (over 60 to date), the majority of which 
do show values of A close to zero, pretty much like the 
planetary bodies orbiting our Sun, although a consider- 
able fraction (nearly 40%) shows substantial misalignment 
( |Albrecht et al.||2012[ ). 

In this paper and in a companion Letter (Paper II of 
the series) by Desidera et al. (2013), we present the first 
results obtained in the framework of the project Global 
Architecture of Planetary Systems (GAPS), a large obser- 
vational programme with HARPS-N which has recently 
been competitively awarded long-term status at the TNG. 
GAPS is a structured, largely synergetic observational pro- 
gramme specifically designed to maximise the scientific re- 
turn in several aspects of exoplanetary astrophysics, taking 
advantage of the unique capabilities provided by HARPS- 
N. The GAPS programme is composed of three main el- 
ements, including a) radial-velocity searches for low-mass 
planets around stars with and without known planets over 
a broad range of properties (mass, metallicity) of the hosts, 
b) characterisation measurements of known transiting sys- 
tems, and c) improved determinations of relevant physical 
parameters (masses, radii, ages) and of the degree of star- 
planet interactions for selected planet hosts. 

In particular, within the framework of the GAPS pro- 
gramme element devoted to investigating the outcome of 
planet-disc and planet-planet interaction scenarios in exo- 
planet systems, we present here the first determination of 
the RM effect for the recently disc overed transiting sys - 
tem Qatar-1 (= GSC 04240-00470) ( |Alsubai et aT1[201l| ). 
This system contains a hot Jupiter orbiting a V=12.84, 
metal-rich K-dwarf star, one of the faintest around which 
a planet has been discovered so far by ground-based sur- 
veys. Furthermore, Qatar-1 represents an interesting study 
case for investigating the star-planet interaction and will 
set constraints on theories of tidal evolution for other sys- 
tems that contain very hot Jupiters orbiting low-mass stars. 
For this purpose, it is mandatory to rely on refined orbital 
parameters as well as on accurate determinations of the 
physical properties for the star and the planet. 

The plan of the paper is the following: in Sec. [2] we 
present the HARPS-N observations and in Sec. [3] we in- 
troduce the complementary photometric observations. In 
Sec. [4] we describe the analysis and present the results of 
the photometric and spectroscopic data, in Sec. [5] we de- 
rive the atmospheric properties of the host star based on 
different methods, while in Sec.[6]we infer the system prop- 
erties related to rotation and activity indicators and discuss 



1 see the RM encyclopedia: 
http://www.aip.de/People/RHeller 
and the TEPCat url: 



http://www.astro.keele.ac.uk/jkt/tepcat/rossiter.html 
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them in terms of star-planet tidal interaction. In Sec. [7] we Table 1. HARPS-N RV measurements of Qatar-1. 
summarise our results and main conclusions. Barycentric Julian Dates (BJD) in TCB units (SI units). 



2. HARPS-N observations and data reduction 

The spectroscopic observations of the transit were obtained 
on 2012 September 3, using the HARPS-N (High Accuracy 
Radial velocity Planet Searcher-North) spectr ograph at 
the 3.58m Telescopio Nazionale Galileo (TNG) (Cosentino 
et al.|2012| ). HARPS-N, a near twin of the HARPS instru- 
ment m operation on the ESO 3.6m telescope at La Silla 
(Chile), covers the wavelength range from 3800 A to 6900 A 
with a resolving power of R^115,000. Each resolution ele- 
ment is sampled by 3.3 CCD pixels. 

While out-of-transit data were obtained with the simul- 
taneous Th-Ar calibration, the in-transit spectra were ac- 
quired with the second fibre on-sky to avoid the risk of 
contaminating the stellar spectrum by the calibration lamp. 
Knowing the spectroscopic orbital parameters, in particu- 
lar the semi-amplitude K of the radial velocity (RV) curve 
and the systemic RV 7 is mandatory for a correct interpre- 
tation of the RM effect. Following the successful acquisition 
of a spectral time series covering the Qatar-lb transit, we 
were prompted to gather additional HARPS-N data aim- 
ing to cover out-of-transit phases and improve the orbital 
solution. 

Thanks to the flexible scheduling of observations inside 
the GAPS programme, during following nights (September 
5, 6, 7, 8, 9, and 11) we were able to obtain seven additional 
spectra evenly distributed over the different orbital phases. 
The RV measurements of Qatar-1 are reported in Table[T] 
All spectra were acquired with an exposure time of 900 s. 
The spectrograph is equipped with an exposure meter to 
accurately measure the flux-weighted mean time of each 
exposure. The RV measurements and corresponding errors 
were obtained using the HARPS-N on-line pipeline, based 
on the numerical cros s-correlation function (CCF) method 
( Baranne et al.||1996[) with the weighted and c/eaned-mask 
modification QPepe et al. 20021, by applying the K5 mask. 

To perform the detailed analysis of light curves and 
radial velocities data the time-tag of each exposure was 
reduced to the sol ar system barycentr ic time using the 
software Tempo2 ( jHobbs et all, |2006|) with the DE405 
JPL Ephemerides ( Standish||1998p in barycentric coordi- 
nate time (TCB) scale (SI units). The assumed celestial co- 
ordinates of the source are: RA(J2000) = 20' 1 13 m 31 s .615 
DEC(J2000) = +65° 09' 43".48, with proper motions (in 
mas/yr) fi a = 7.1, /is — 58.0. In Tablefll we report the RV 
after barycentering with Tempo2 in TCB units. 

3. Ancillary data: transit photometry 

New photometric data of Qatar-lb transits were obtained 
with the Asiago 1.82m and Calar Alto 1.23m telescopes. 
The journal of photometric observations is given in Table[2] 
while the photometric data are provided in Table|3] 



3.1. Photometric observations from the Asiago Observatory 

Two complete transits of Qatar-lb were observed on 2011 
May 29 and 2012 Aug 24 within the TASTE project 



BJD (TCB) 


RV 

[ms- 1 ] 


error 

[ms- 1 ] 


S/N t 


t 


2456174.4128054 


-38006.8 


4.7 


20 


b 


2456174.4247685 


-38011.8 


4.4 


21 


b 


2456174.4355706 


-38021.1 


4.4 


21 




2456174.4463817 


-38016.3 


4.4 


22 




2456174.4572017 


-38049.8 


4.7 


20 




2456174.4680268 


-38061.4 


4.7 


20 




2456174.4788509 


-38093.0 


4.5 


21 




2456174.4896710 


-38102.7 


4.6 


20 




2456174.5005001 


-38097.8 


4.8 


20 




2456174.5113162 


-38115.5 


4.8 


20 


b 


2456174.5232873 


-38128.7 


5.0 


19 


b 


2456176.5702612 


-38084.6 


7.5 


13 





2456177.5788490 


-38311.6 


7.4 


14 





2456177.5954171 


-38310.1 


10.4 


11 





2456178.5619923 


-37885.3 


13.5 


9 





2456179.6162127 


-37836.5 


19.7 


7 





2456180.5669062 


-38291.1 


6.7 


16 





2456182.5996329 


-37781.7 


6.0 


17 






f: signal-to-noise ratio per pixel at 5500 A; 
f: i = in-transit, = out-of-transit, b = data taken just 
before or after the transit, which are used for both the 
orbital and the RM fit. 



the autoguide. Both observations were performed with the 
Asiago Faint Object Spectrograph and Camera (AFOSC) 
at the 1.82m Copernico telescope in northern Italy. AFOSC 
is a classical focal-reducer camera equipped with a thinned, 
back-illuminated E2V 42-20 2048 x 2048 CCD (0.26"/pix 
in unbinned mode) , providing a 9' x 9' field of view. Both 
transits were observed in imaging mode through a standard 
Cousins R filter, with a constant exposure time of 7 s. CCD 
windowing and 4x4 binning were set to increase the time 
series duty-cycle, while a suitable set of reference stars was 
imaged on the same read-out window. Stellar images were 
defocused to ~ 4" FWHM to avoid saturation and minimize 
flat-field residual errors. 

After a standard correction for bias and flat-field, the 
frames were reduced with the STARSKY code the TASTE 
photometric pipeline (Nascimbeni et al. 2013), which also 
includes a diagnostic tool to discard the data points most 
affected by transparency variations. Differential light curves 
were extracted by normalising the raw flux from the target 
with an optimally weighted average flux from the reference 
stars. Finally, light curves were corrected for systematic er- 
rors by decorrelating flux against external parameters (in- 
cluding e.g., star position on the detector, FWHM, back- 
ground level, and airmass) and selecting the solution with 
the smallest off-transit scatter. 



(Nascimbeni et al. 2011). The weather conditions in both 



nights were characterised by veils and thin cirrus, yet nei- 
ther caused interruptions in the time series nor stops in 



3.2. Photometric observations from the Calar Alto 
Observatory 

Three transits of Qatar-lb were observed on 2011 August 
25 and on 2012 July 21 and September 10, using the 1.23 m 
telescope at the German-Spanish Calar Alto Observatory 
(CAHA) near Almeria (Spain), which was already success- 
fully used t o follow-up several transiting planets (Mancini 
et al. 2012). During the 2011 observations, we used the 



Covino et al.: GAPS. I. RM effect in Qatar- 1 

Table 2. Details of the photometric observations presented in this work. iV b s is the number of observations, Moon is 
the fractional illumination of the Moon at the midpoint of the transit, and i e xp is the exposure time in seconds. The 
aperture sizes are the radii in pixels of the software apertures for the star, inner sky, and outer sky, respectively. Scatter 
is the r.m.s. scatter of the data versus a fitted model in mmag. Times and dates are in UT. 



Obs. 



Camcra/CCIT 



Date 



UT(Start->End) JV P 



Filter 



Airmass Moon Aperture Scatter 



Asiago 
CAHA 
CAHA 
Asiago 
CAHA 


AFOSC/E2V 42-20 
SITE#2b/2kx2k 
DLRMKIII/4kx4k 
AFOSC/E2V 42-20 
DLRMKIII/4kx4k 


2011 05 29 

2011 08 25 

2012 07 21 
2012 08 24 
2012 09 10 


22:54 -» 02:19 
23:31 -* 04:02 
20:29 -> 00:32 
23:09^ 02:13 
00:20 -» 03:41 


1218 
186 
112 

1134 

78 


7 

60 

120 

7 
120 


Cousins R 
Cousins R 
Cousins R 
Cousins R 
Cousins R 


1.11 -> 1.06 
1.14 -» 1.78 
1.45 -> 1.13 
1.09 -> 1.38 
1.30 — ^ 2.02 


7% 
61% 
77% 
55% 
77% 


10,12,22 
12,30,45 
10, 30, 50 

5,6,16 
21,33,50 


1.25 a 
1.65 
0.78 

1.41 a 
0.87 


< a > After 


binning 





















Table 3. Excerpts of the R-band light curves (LC) of 
Qatar- 1: this table will be made available at the CDS. A 
portion is shown here for guidance regarding its form and 
content. 



Tclcscopc/LC BJD (TCB) 



Diff. mag. Uncertainty 



Asiago #1 
Asiago #1 

Asiago #2 
Asiago #2 

CA 1.23m #1 
CA 1.23m #1 

CA 1.23m #2 
CA 1.23m #2 

CA 1.23m #3 
CA 1.23m #3 



2455711.4545895 
2455711.4548095 

2456164.4647472 
2456164.4648472 

2455799.4805145 
2455799.4814715 

2456130.3541082 
2456130.3561482 

2456181.5142293 
2456181.5159694 



0.0016630 
0.0052700 

0.0025770 
0.0021340 

-0.0002695 
0.0002191 

-0.0004088 
0.0007516 

-0.0000350 
-0.0001722 



0.0034607 
0.0034568 

0.0027072 
0.0027099 

0.0014395 
0.0014412 

0.0008100 
0.0007696 

0.0005951 
0.0005738 



2kx2k SITE#2b optical CCEQ with a FOV of 16' x 16' 
and a pixel size of 24 p,m, which translates to a pixel scale 
of 0.5" per pixel. We defocused the telescope in order to 
lower the flat-fielding noise. The photometric data were 
gathered through a Cousins R filter with an observing ca- 
dence of 60 sec. To limit the dead time between exposures, 
we reduced the amount of time lost to CCD readout by 
reading out only a small window. Autoguiding was used. 
The two 2012 transits were obtained in filter Cousins R 
with the new DLRMKIII camera, which is equipped with 
an E2V CCD231-84-NIMO-BI-DD sensor with 4kx4k pix- 
els of 15/xm and a FOV of 21' x 21'. Both transits were 
observed with the CCD unbinned, the telescope heavily de- 
focused, and an observing cadence of 120 s. The autoguider 
was properly focused to preserve a good pointing of the 
telescope. 

T he observations we re analysed using the idl pipeline 
from iSouthworth et al.| (|2009t) . The images were debiased 

methods, then subjected to 



( [20091) . 
mdard 



Southworth et al 
and flat-fielded using stancH 

aperture photometry using the IDL task APEr|^J Pointing 
variations were followed by cross-correlating each image 
against a reference image. We chose the aperture sizes and 
comparison stars that yielded the lowest scatter in the fi- 
nal differential-photometry light curve. The relative weights 
of the comparison stars were optimised simultaneously by 
fitting a second-order polynomial to the outside-transit ob- 
servations to normalise them to unit flux. 



2 This CCD has since been decomissioned. 

3 Part of the "ASTROLIB subroutine library" 



distributed by 



4. Data analysis 

In this section we describe the analysis of the data pre- 
sented in Sec. [3] and present the results we obtained. The 
analysis of our five photometric data sets was performed by 
employing the software code JKT EBOP (version 28 ) to fit 
a transit light-curve (LC) model ( Southworth|2008 1 . 



4.1. New ephemerides for Qatar- 1 b transits 

As a first step we derived new ephemerides by combining 
the determinations of the mid-transit times from our five 



data sets with those derived from |Alsubai et aL (2011 1. All 
timings were placed on BJD (TCB) time system. The re- 
sulting measurements of transit midpoints were fitted with 
a straight line to obtain new orbital ephemerides: 



T = 2 455 518.41094(28) + 1.42002504(71) x E, 



(1) 



where E is the number of orbital cycles after the refer- 
ence epoch (the mi dpoint of the first transit observed by 
Alsubai et al. 2011) and the quantities in brackets denote 
the uncertainty in the last digits of the preceding number. 
The corresponding O— C diagramme is shown in Fig. [II in 
which the mid-transit times available from the Exoplanet 
Transit Database (ETD Pj are also displayed, though they 
were not used in the fit. 



4.2. Combined solution of photometric light-curves 

The main parameters of the model are the fractional radii of 
the star and planet, r* and r p , defined as the stellar radius 
R+ and the planetary radius R p scaled by the semi-major 
axis a, respectively, and the orbital inclination i p . 

The LC solution was attempted using the linear and 
quadratic limb-darkening (LD) law, and with the LD coef- 
ficients either fixed at the theoretical values tabulated by 
Claret & Bloemen (2011) or included as fit parameters. 
We obtained the best result by fitting the LD(ui) coeffi- 
cient using a linear law. We fitted individual light curves 



for the sum and ratio of the fractional radii r p + r+, r p /r+, 
the orbital inclination i, the central transit time To, and 
the limb darkening coefficient LD(ui). Evaluation of robust 
uncertainties on our best-fit parameters was performed by 
a bootstrap algorithm, i.e. generating for each light curve 
10 000 resampled data sets to be fitted, and an alysing the 
resulting distribution (Southworth et al. 2005). The light- 
curves and their best-fitting models are shown in Fig.[2[ 
whereas the results of our analysis are reported in TableH] 



NASA: idlastro.gsfc.nasa.gov/ 



http://var2.astro.cz/ETD/ 
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oASai a AEe 

Fig. 1. O— C diagrammc for Qatar-1 transit times. The small symbols in colour refer to the Asiago (green), Calar Alto 
(blue) and Alsubai et al. (2011) (red) data used for redetermination of the ephemerides. Open circles represent transit 
times available from the ETD. 



along with the weighted means (WMs) for each parameter 
evaluated over all five transits. For some parameters (r p /r+ 
and LD(ui) the best-fit values are more scattered than ex- 
pected by applying Gaussian statistics. A possible reason 
for these slight discrepancies could be intrinsic variations 
due to stellar activity of Qatar-1. Stellar spots, faculae and 
other active regions are known to alter the parameters in- 
ferred from photometric data, even when they are not oc- 
culted during the transit and thus no unusua l feature is 
visible in the light curves (Bal lerini et al.|[2012 ). The ratio 
r p /r+ is one of the most critical parameters in this regard. 
We checked for long-term variability of Qatar-1 by perform- 
ing differential photometry with the same comparison stars 
on the Asiago time series, Al and A2, taken at the same 
site and airmass. We found a brightness variation between 
the two epochs of +0.034 ± 0.005 and +0.030 ± 0.006 mag 
(using two different reference stars), while the variation of 
our "check" star is -0.004 + 0.004 mag. Although this find- 
ing is not fully conclusive, it suggests that stellar activity 
could play a role and possibly explain the differences in the 
LC solutions for different epochs. 
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4.3. Improved spectroscopic orbit solution 

The 11 RV measurements corresponding to out-of-transit 
phases (see Table[l]) form the data set we used to fit the 
spectroscopic orbit solution. 

The solution of the RV curve was performed considering 
both the cases of circular and eccentric orbit. The results 
of our RV curve fits are summarised in Table [5J and the 
orbital solution corresponding to the fit obtained with the 
eccentricity as free parameter is shown in Fig. [3] The er- 
rors on the best-fit parameters were determined with the 
bootstrap method from 1000 mock data sets. 

The comparison with the best-fit parameters reported 



Alsubai et al. (2011) shows a number of significant differ- 



ences. The Alsubai estimate of 7=-37835±63 ms^ 1 deviates 
from ours formally by more than 3c. The most likely ex- 
planation for this discrepancy lies in a systematic RV zero- 
point difference between the two instruments used. The 
higher precision of the HARPS-N RV measurements allows 
us to rule out an eccentric orbit for Qatar-lb. Our deter- 
mination of e = 0. 02010°^ is compatible with a circular 
orbit within 2a, as expected for clos e-in planets with or - 
bital periods shorter than a few days (Husnoo et al.|2012|. 



Fig. 2. Photometric transit light-curves of Qatar-1 and cor- 
responding residuals between data points and best-fit solu- 
tion. 



Therefore, for the following analysis we adopted the best-fit 
parameters from the circular solution. 

Finally, we obtained K=265.7 ± 3.5 ms" 1 , which is 



about 20% higher than the value in Alsubai et al. (2011). 
Accordingly, the estimates of the mass and density of 
Qatar-lb will be higher by the same percentage. 

4.4. Analysis of the RM effect and determination of the 
spin-orbit alignment 

To analyse the RM effect we implemented a numerical 
model based on the following assumptions. We tried to re- 
produce the observed CCF by modelling the average photo- 
spheric line profile. The stellar disc is sampled by a matrix 
of 2000x2000 elements, each element being represented by 
a Gaussian line profile with a given width u e i, Doppler- 
shifted according to the stellar rotation, and weighted by 
appropriate limb-darkening coefficients. The resulting line 
profile is then convolved by the instrumental profile (as- 
sumed to be Gaussian) of HARPS-N, oip = 1.108 km s _1 . 
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Fig. 3. HARPS-N radial velocity data and best-fit spectro- 
scopic orbit solution for Qatar-1. Upper panel: RV curve 
versus BJD showing only the out-of-transit data used for 
the orbital solution (solid red line), and the correspond- 
ing residuals. Lower panel: phase-folded RV curve (solid 
red line) and data used for the orbital fit drawn as black 
symbols. The RV measurements corresponding to in-transit 
phases are also displayed as light-blue symbols. 




— D 

1 o 
o 
6 -D 

-D 



H 1 1 1 



I i * i t f ; j 



6174.42 6174,44 6174.46 6174.48 6174.50 6174.52 

BJD-2450000 



Fig. 4. RM effect best-fit solution (solid, red line) with 
RV measurements (black star symbols) and computed RV 
values (red open circles), yielding the O— C beneath. The 
Keplerian radial velocity curve is also drawn as a blue dot- 
ted line. 



Table 5. Orbital parameters derived from the spectro- 
scopic orbit solution of Qatar-1. The parameters that were 
kept fixed in the fit are in italics. 

Orbital parameter 

Eccentric fit Circular fit 
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The model also takes into account the actual occulted area 
of the stellar photospheric disc and the smearing due to the 
planet's displacement during an exposure. The correspond- 
ing RV shift is then computed by a Gaussian fit of this re- 
sulting line profile, analogously to the HARPS-N pipeline. 
The model depends on twelve parameters. The orbital pe- 
riod Porfc, mid-transit epoch To, transit duration T14, star 
and planet radii R* and R p , the linear limb-darkening co- 
efficient LD(ui), and the impact parameter b were fixed to 
the values we derived from the light-curve solution, while 
the star RV semi-amplitude K was set to the best-fit value 
obtained from the spectroscopic orbit with null eccentric- 
ity. The remaining parameters, i.e. the stellar projected ro- 
tational velocity vsin J, the systemic velocity 7, the stellar 
disc resolution element line width a e i , and the orbital obliq- 
uity A were treated as free parameters. 

By usi ng a "trust region" least-squares mini misation 



algorithm ( Branch et al.||1999 
tained the following best-fit v 
vsin/ = 1.7 ± 0.3kms _1 , a e i 



Byrd 

aiues: 



et al. 



A 



1988) 



we ob- 
7.1 deg, 

"\7 = 
The er- 



2238 ± 155 ms- 
-38059.5 ± 2.0ms" 1 , and normalised x 2 = 1-217. 
rors were determined with the same bootstrap method as 
for the orbital RV curve fit from 200 mock data sets. 

Although an independent determination of the systemic 
velocity 7 is provided by the orbital RV curve fit, we pre- 
ferred to treat 7 as a free parameter, because of its strong 
correlation with the parameter A (see Fig. [5]). In fact, by 
fixing the value of 7 to the value derived from the orbital 
RV curve fit, we obtained a best-fit value of A = —1°, with a 
normalised \ 2 of 1.301. By letting 7 free to vary, we also al- 
lowed for a possible RV shift caused by stellar spots, whose 
configuration on the stellar disc can be assumed to remain 
constant over the time-span of the transit. 

This analysis shows that the system is well aligned. The 
results are reported in Tablerol Fig. ^ shows the x 2 -maps for 
various pairs of parameters. 

A preliminary calibration of the FWHM of the CCF 
in terms of projected rotational velocity based on a sam- 
ple of stars with planets observed as part of GAPS yields 
vsin/=1.6±0.5kms _1 , consistent with the value derived as 
part of the RM effect modelling. 



5. Characterisation of the host star 

A proper determination of the stellar parameters is manda- 
tory to obtain the physical properties of the planet. 
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Fig. 5. x 2 -maps showing the correlation between different couples of parameters. 



Because of the uncertainty on the amount of interstellar 
extinction, photometric data alone cannot provide an ac- 
curate determination of the fundamental paramete rs of the 
star. Accord ing to the TASS Mark IV catalogue (Droege 



et al 



20061), the apparent V-band magnitude of Qatar - 



84 mag. The K2V spectral type of the host 
star (from Alsubai et al. 20111 translates into an abso - 
lute magnitude of M„=6.5 mag ( |Straizys fe Kuriliene 1981 ). 
Depending on the amount of reddening, the distance to 
Qatar-1 is expected to span roughly between ~190pc for 
negligible extinction, and ~130 pc for a normal interstellar 
extinction law (i.e., R v = 3.1) and the maximum colour 
excess value (i.e., Eb-v = 0.232 ±0.003 ) obtained from the 



reddening map by Schlegel et al. (1998). 



5.1. Spectroscopic determination of stellar parameters 

The HARPS-N spectra were used to perfom a spectroscopic 
characterisation of the host star Qatar-1 and estimate the 
effective temperature T c ff, surface gravity logg, projected 
rotational velocity vsin J, and the iron abundance [Fe/H] ) . 
The effective temperature was initially determined from 
the HARPS-N spectra by applying the method of equiv- 
alent width (EW) ratios of spectral absorption lines by 
means of the ARES^] automatic code ( Sousa et al. 2007 ), 
using the calibration for FGK dwarf stars by Sousa et al. 
(2010). For this purpose, we used an average spectrum 
of Qatar-1 obtained by a weighted mean of the 18 spec- 
tra available (after verifying that no contamination was 
present in the spectra taken with the simultaneous ThAr 



calibration), each properly shifted by the corresponding 
radial velocity to the rest wavelength frame, obtaining 
T cff = 4990 ± 100K. 

Furthermore, the atmospheric s tellar paramet ers were 
derived using the program MOOG (Snedcn 1973), version 
2010, and through EW measurements of iron lines as de- 
scribed in detail by Biazzo et al. (2012). In particular, the 



effective temperature was determined by imposing the con- 
dition that the Fei abundance does not depend on the ex- 
citation potential of the lines, the microturbulcncc velocity 
(£) was derived by imposing that the surface Fe I abundance 
is independent on the line EWs, and the surface gravity 
was estimated by imposing the Fe i/Fe n ionization equilib- 
rium. Then, we were also able to measure the iron abun- 
dance of Qatar-1. This analysis was performed differentially 
with respect to the Sun. For this purpose, we analysed the 
Qatar-1 spectrum and a Ganymede spe ctrum acquired with 
HARPS at ESO ( [Mayor et aT]|2003[ ). We thus obtained 
logn(Fei)© = logn(Feli) = 7.53±0.05 for the Sun. In 
the end, we found the following atmospheric parameters 
and iron abundance for the star: T e s = 4820±100 K, logg = 
4.43 ±0.10, £ = 0.90 ±0.05 kms" 1 , [Fei/H]= 0.15 ± 0.10, 
and [Fen/H]= 0.1 5 ± 0.08, which agre e with the previous 
determinations by Alsubai et al. (2011). 



http://www.astro.up.pt/~sousasag/ares/ 



In addition, the stellar parameters were derived inde- 
pendently using the methodo logy presented in |Santos et al.| 
(2004) and Sousa et al. ( |2008 ). This analysis was also based 
on iron line excitation and ionization equilibrium, using a 
grid of Kurucz (1997) model atmospheres and the 2002 ver- 
sion of MOOG. The ARES code was used to measure the 
equivalent widths for a sub-set of iron lines from the list in 
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Fig. 6. Portion of Qatar-1 average spectrum (black solid 
line) and synthetic spectrum (red solid line) couplcorre- 
sponding to the following set of (fixed) parameters: T c g — 
4910K, \ogg = 4.66, £ = OJOkms" 1 , and wsin/ = 
1.7 km s -1 . The positions of five iron lines and two lines 
of titanium and calcium are marked. The inset shows a 
zoom on the spectral window containing three Fe I lines, 
and the comparison with synthetic spectra corresponding 
to [Fei/H]=0.10 (green, long-dashed line), 0.20 (red, solid 
line), and 0.30 (blue, dotted line). 



Sousa et al. (2008) that is especially suited for determin- 
ing the parameters of stars with temperature below 5200 K 
(for details see Tsantaki et al., A&A, submitted). The 
derived parameters are T of f=4786±95K, log#=4.41±0.24, 
£ = 0.78 ± 0.29 kms" 1 , and [Fe/H]=0.18±0.06. These val- 
ues agree fairly well with those reported above. 

As a final consistency check, we exploited the fact that 
there is only one temperature for which the surface grav- 
ity from evolutionary models will meet the one obtained 
from the ionization equilibrium. By imposing this con- 
dition, the surface gravity was obtained by comparison 
with the BaSTJj evolutionary models. This yields the set 
of values T cS = 4910 ± 100 K, logg = 4.66 ± 0.10, and 
[Fe/H]= 0.20 ± 0.10. Fig.|6] shows the observed and the 
synthetic spectrum corresponding to the latter set of pa- 
rameters and adopting for the projected rotational veloc- 
ity the value vsinl = 1.7kms _1 , as derived from the RM 



best-fit solution (see Sect 4.4 ) . This set of values was finally 



adopted by us in the following. 



5.2. SED, reddening, and distance 



We constructed the spectral energy distribution (SED) 
of Qatar-1 by merging the V and R optical magnitudes 



from the TASS Mark IV catalogue (Droege et al. 2006) 



with the i nfrared photometry from the 2MASS and WISE 



databases ( |Cutri et al.|2003 Cutri fe et al.|2012 ) , as shown 
in Fig. [7] We simultaneously fitted the colours encompassed 
by the SED to ad hoc synthetic ma gnitudes derived from a 
NextGen stellar atmosphere model (Hauschil dt et al.1J 999) 
with the same T g, logg, and metallicity as the star. To 
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Fig. 7. De-reddened spectral energy distribution of Qatar- 
1, for A v — 0.1 and d=195pc. Optical V and R pho- 
tometr ic data are taken f rom the TASS Mark IV cata- 
logue (|Droege et al. 2006). Infrar ed data are taken from 



the 2MASS a nd WjffEliatabases ( [Cutri et al.|[2^0^[Cutri 
I fe et al.|2012 ). The NextGen synthetic low-resolution spec- 
trum ( |Hauschildt et al. 1999) with the same photospheric 
parameters as Qatar-1 is over-plotted with a light-blue line. 
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Fig. 8. Portion of the HARPS-N averaged spectrum of 
Qatar-1 highlighting the weak interstellar Nai lines. 



estimate the interstellar extinction A v and distance d to 
Qatar -1, we followed the method described in |Gandom et al.| 
( 2008 ) and combined the availabe set of photometric data 
with the spectroscopically derived parameters. Assuming a 
total-to-selective extinction R v = A v /Eb-v = 3.1 and a 
black body emission at the star effective temperature and 
radius yields an extinction A v — 0.10 ±0.10 mag and a dis- 
tance to the star d = 195 ± 25 pc, fairly consistent with the 
range of values estimated above. 

Weak interstellar components are visible within the 
broad and blended Nai Di j2 stellar doublet, as shown 
in Figure[8J Their equivalent widths are 0.017 and 0.009 
(±0.001) A^ respectively, i.e. very close to the asymptotic 
2.0 limit ratio for optically thin 5889.951 and 5895.924 A 
lines. Some contamination from the adjacent not perfectly 
subtracted night -sky lines is possible. Adopting the Munari 
& Zwitter ( 1997 ) calibration between reddening and equiva- 



URL |http://wwwas.oats.inaf.it/IA2/BaSTI/ 



lent width of interstellar Nai 5889.951 A line, the reddening 
affecting Qatar-1 is £:b-v=0.004 (±0.0005). 
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Finally, using a method based on the near-infrared 
SED publishe d in Masana et al. (20061, as described in 
Section 2.3 of Ribas et al. (2010), and adopting the values 



of logg and [Fe/HJ derived from the spectroscopic analy- 
sis (though their impact on the final values is very small, 
similarly to the immunity of the infrared flux method to 
those parameters), we obtained an effective temperature of 
T e ff = 4880 ± 70 K. The main source of error is the uncer- 
tainty in the extinction value, assumed to be Av=0.1+/-0.1 
mag. Hence, the resulting effective temperature perfectly 
agrees with the spectroscopic determination. 

5.3. Star and planet physiscal properties 

Table[6] summarises the physical properties we derived for 
the star and the planet in the Qatar- 1 system. To de- 
termine the stellar fundamental properties, we exploited 
quantities obtained directly from the analysis of the tran- 
sit ligh t curves following the prescriptions by Soz zetti et al.| 
( |2007 ). In particular, we derived the density of the star p+ = 
1.62 ± O.O8y00 from the scaled stellar radius R+/a and used 
it together with T off and [Fe/H] obtained from the spec- 
troscopic anaslysis to infer the stellar mass and radius by 
comparison with stellar evolution models. Using the BaSTI 
models, this yields a stellar mass M* = 0.85 M ± 0.03 M Q 
and a radius i?* = 0.80 ± 0.05 -R ^jwhi ch coincide with the 



values reported by Alsubai et al. ( 2011 ). In turn, the surface 
gravity of the planet, logg p = 3.372 ±0.024, and the planet 
equilibrium temperature, Tp — 1389 ± 39 K, with their cor- 
responding uncer tainties were obtained following |Sozzetti 
et al.| ( |2007[ ) and |Cowan fc Agol| poTTj ), respectively. The 
planetary radius and mass are found to be 1.18 ± 0.09 Rj 
and 1.33 ± 0.05 Mj. Hence, while only marginally larger in 
radius, the planet is found to be signifi cantly more massive 
than reported by Alsubai et al. (2011). 



6. Discussion 



As discussed by Enoch et al. (2012), the radius of an exo- 
planet may be affected by a number of factors, which in- 
clude the mass of the planet and its heavy element content 
as well as the irradiation received from the host star. For 
Qatar-1 b, we expect strong irradiation and tidal effects be- 
cause of its proximity to the host star. Below, we examine 
rotation and activity indicators and discuss the possible 
evolutionary effects caused by tidal interaction. 

6.1. Stellar rotation, activity, and age 

Adopting the vsinl of 1.7 ± 0.3kms _1 and radius of 
0.80 ± 0.05 R© for Qatar-1, and assuming that the star is 
viewed equator-on, the rotation period is calculated to be 
24 ± 6 days. The gyrochronological a ge, t avro , can h ence be 



estimated using the relation given by Barnes ( 2007 ) (equa- 
tion 3). For a (B — V) = 0.9, corresponding to a T c g of 
4910 K of the star, we obtain t gyro = 1.7 ± 1.1 Gyr. 

On the other hand, adopting the revised gyrochronology 
by |Lanza| ( ]2010[ ) for stars with hot Jupiters, we estimate an 
age of ~ 4.5 Gyr that better agrees with that derived from 
stellar evolutionary tracks (cf. Fig. 5 of Alsubai et al. 2011). 

Chromospheric emission is present in the Can H&K 
lines. Unfortunately, the S/N of individual HARPS-N spec- 
tra in the wavelength range below 4000A is inadequate 



to obtain reliable measurements of the log R' H K index. 
Therefore it is not possible to look for eventual variations 
in the chromo spheric activity . Howe ver, using the passband 
definitions in Duncan et al. ( 1991 ) and a preliminary cal- 
ibration of the flux in the Call H&K lines, we obtained 



Shk = 0.389 from the average spectrum used in Sec. |5.1 
Then, by adopting (B — V) = 0.9 and the calibrations of 
Noyes et al.] (l984[), we find logR^ K = -4.60. Therefore 



Qatar-1 appears to be a mod erately active star, si milar to, 
e.g., HD 189733 and TrES-3 dKnutson et aI1|2010l). 

Using the relations g i ven b y |Noyes et al. (1984) and 
Mamajek & Hillenbrand (2008), the expected rotation pe- 
riod for logR^ = -4.60 and (B - V) = 0.9 would 



be 23.8 d and 23.0 d, respectively, hence within the range 
estimated from vsinJ. However , using relation (3) from 
[Mamajek & Hillenbrand (20081 of age as a function of 
logRy K , yields an age of about 1.1 Gyr, which is consider- 
ably you nger than t he value inferred from gyrochronology 
(see also Pace 2013 for an updated analysis on the use of 
chromospheric activity as an age indicator). 

Finally, we note that Qatar-1 b fits well with i n the 



log R^ K -log g p correlation pointed out by [Hartman (2010) 
for close-in planets (a < 0.1 AU) with M p > 0.1 Mj, or- 
biting stars with 4200 K< T cff < 6200 K. It ap pears also 
to be con s istent with the indication reported by |Krejcova 
& Budaj ( J2012 ) that the level of the chromospheric ac- 
tivity of stars with T e g < 5500 K hosting close-in planets 
(a < 0.15 AU) may be enhanced by the presence of Jupiter- 
mass planets. 



6.2. Star-planet tidal interaction 

The orbital period of the planet in the Qatar-1 system 
is much shorter than the rotation period of the star es- 
timated in Sec. 6.1 Therefore, tides produce a decay of the 
orbit with a continuous transfer of angular momentum from 
the orbital motion to the spin of the star. The orbital an- 
gular momentum is only 0.28 of the present stellar spin, 
which is insufficient to reach an equilibrium synchronous 
state. In other words, as a consequence of the tidal evo- 
lution, the planet is expected to be engulfed by the star. 
The timescale for the orbital decay t^ 1 = (l/a)(da/dt), 
where a is the orbital semi-major axis and t the time, de- 
pends on the efficiency of the tidal dissipation inside the 
star that is parameterised by its modified tidal quality fac- 
tor Q' s . Considering the tidal dissipation efficiency required 
to account for orbital circularisation and sy nchronisation in 
close binary systems, lOgilvie & Lin (2007) estimated that 
Q' s is of the order of 10° for late-type mam-sequence stars. 
Adopting this value and the stellar and planetary param- 
eters in Table|6j we obtain r a ~ 0.5 Gyr for Qatar-1 using 



the tidal evolutionary model of Leconte et al. (20101. The 
timescale for the alignment of the orbital angular momen- 
tum and the stellar spin can be estimated from the same 
tidal model and is r c _1 = (l/e)(de/dt), where e is the obliq- 
uity of the stellar equator with respect to the orbital plane. 
Assuming Q' s = 10 6 and an initial obliquity e(t — 0) = 30° 
yields r e rj 0.2 Gyr. 

It is interesting to note that the timescale of orbital de- 
cay computed with Q' s = 10 6 is significantly shorter than 
the estimated age of the star from Sect |6.1| The difference 
is not as dramatic as in the case of, e.g., OGLE-TR-56 



( |OgilviefcLm|2007[ ), or Kepler- 17 ( |Bonomo fc Lanza|2012| ), 
where the orbital decay timescale turns out to be as short as 
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40-70 Myr, but this result suggests that the use of the same 
value of Q' s as derived for close stellar binary systems also 
for star-planet systems is not appropri ate (cf. other anal- 
yses trying to constrain Q' s , e.g., |Patzold fc Rauer 2002 
Ogilvie fe Lin|20071 [Carone fc Patzold||2007|| Jackson et al. 



2009 Patzold et al.|2012 ) . In principle, a lower limit on Q 
can be established by an accurate timing of the transits ex- 
tended over a baseline of several decades, which allows one 
to test the predictions of the above models. Specifically, for 
Q' s = 10 6 , the orbital decay is expected to produce a varia- 
tion of the observed epoch of mid-transit with respect to a 
constant-period ephemeris of O — C ~ 17 s in twenty years. 
The tidal evolution o f stellar obliquity has been recently 
discuss ed by Lai (|2012) from a theoretical point of view, 
and by Albrecht et al. (20121 considering RM observation 



statistics. 



7. Summary and conclusions 

We reported on observations of the RM effect for the Qatar- 
1 system and a new determination of the orbital solution 
based on HARPS-N at TNG observations. With these data 
we also obtained the spectroscopic characterisation of the 
host star. Combining the radial velocity data with new tran- 
sit photometry and with photometric data from the litera- 
ture, we derived new ephemerides and improved the orbital 
parameters for the system. The most notable results from 
this analysis can be summarised as follows: 

1. The new spectroscopic orbital solution is found to be 
consistent with a circular orbit and any significant ec- 
centricity of the orbit can be ruled out. 

2. The RM effect was measured and the sky-projected 
obliquity was determined, from which we can conclude 
that the orbital plane of the system is well-aligned with 
the spin axis of the star. This result and the derived 
properties of Qatar-lb appear to be in line with the 
general trend observed for close-in pl anets around stars 
cooler than about 6250 K ( |Albrecht et al.||2012[ ). 

3. The host star is confirmed to be a slowly rotating, 
metal-rich K-dwarf, which yet is found to be moder- 
ately active, as inferred from the strength of the chro- 
mospheric emission in the Call H&K line cores and 
from changes in the photometric light-curves at differ- 
ent epochs, which hint at the presence of stellar spots. 

4. The planet is found to be s ignificantly more ma ssive 
than previously estimated by Alsubai et al. (2011). 

5. Qatar- 1 appears to be consistent with the indication 
that the level of the chromospheric activity in stars 
cooler than 5500 K that host close-in giant planets may 
be enhanced. An attempt at estimating the timescale 
for the orbital decay by tidal dissipation was also pre- 
sented, which deserves further investigation. 
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Table 4. Results from the analysis of individual and combined tranist light-curves 





Date/ID 


r*±r p 


r p / r* 


i P 


r* 


r p 


b 


LD (mi) 


2011.05.29/A1 


0.1850±0.0120 


0.1441±0.0041 


83.67±1.18 


0.1617±0.0100 


0.0233±0.0021 


0.682±0.038 


0.49±0.16 


2012. 08. 24/ A2 


0.1716±0.0093 


0.1431±0.0022 


84.99±0.74 


0.1501±0.0077 


0.0215±0.0014 


0.582±0.047 


0.50±0.09 


2011.08.25/C1 


0.1793±0.0089 


0.1465±0.0037 


84.56±0.87 


0.1564±0.0075 


0.0229±0.0016 


0.606±0.068 


0.67±0.12 


2012.07.21/C2 


0.1882±0.0043 


0.1535±0.0009 


83.36±0.35 


0.1632±0.0036 


0.0250±0.0006 


0.709±0.022 


0.37±0.11 


2012.09.10/C3 


0.1834±0.0062 


0.1496±0.0022 


84.04±0.53 


0.1595±0.0051 


0.0239±0.0011 


0.651±0.038 


0.64±0.09 


WM 


0.1841±0.0030 


0.1513±0.0008 


83.82±0.25 


0.1601±0.0025 


0.0241±0.0005 


0.675±0.016 


0.54±0.05 



N.B. 



Rp/a; 



R*/a 



Table 6. Summary of planet and star parameters for the Qatar-1 system. 



Parameter [Units] 



Symbol Value 



Method 



Transit epoch (BJD_TCB-2450000.0) [days] 


To 


5518.41094 ±0.00016 


(1) 


Orbital period [days] 


P 


1.42002504 ±0.00000071 


(1) 


Transit duration [days] 


T u 


0.0678 ±0.0010 


(2) 


Planet/star area ratio 


(Rp/R*) 2 


0.02280 ±0.00011 


(2) 


Planet/star radii ratio 


Rp/R* 


0.1513 ±0.0008 


(2) 


Orbital inclination [degrees] 


i P 


83.82 ±0.25 


(2) 


Impact parameter 


b 


0.675 ±0.016 


(2) 


Scaled stellar radius 


R*/a 


0.1601 ±0.0025 


(2) 


Scaled planet radius 


Rp/a 


0.0241 ± 0.0005 


(2) 


Star reflex velocity [km s _1 ] 


K* 


0.2657 ±0.0035 


(3) 


Orbital eccentricity 


e 





(3) 


Systemic velocity [km s _1 ] 


7 


-38.0558 ±0.0002 


(3) 


Star density [p©] 


P* 


1.62 ±0.08 


(4) 


Planet surface gravity [cgs] 


logffp 


3.372 ±0.024 


(4) 


Star mass [Mq] 


M* 


0.85 ±0.03 


(5) 


Planet mass [Mj] 


M p 


1.33 ±0.05 


(6) 


Orbital semi-major axis [AU] 


a 


0.02343 ±0.0012 


(7) 


Star radius [R©] 


R* 


0.80 ±0.05 


(7) 


Star surface gravity 


logs* 


4.55 ±0.10 


(7) 


Planet radius [Rj] 


Rp 


1.18 ±0.09 


(7) 


Planet density [pj] 


Pv 


0.80 ±0.20 


(7) 


Spin-orbit misalignment [degrees] 


A 


-8.4 ±7.1 


(8) 


Star projected rotational velocity [km s _1 ] 


vsin I 


1.7 ±0.3 


(8) 


Star effective temperature [K] 


T e B 


4910 ± 100 


(9) 


Stellar metallicity 


[Fe/H] 


0.20 ±0.10 


(9) 


Star age [Gyr] 


Age* 


Ri4.5 


(10) 


Planet equilibrium temperature [K] 


T P 


1389 ± 39 


(11) 



(1) Derived by combining the determinations of the mid-transit times from our photometric data sets with those derived from 



Alsubai et al. (2011 I. Kept fixed in the orbital and Rossiter RV fitting. (2) Derived as weighted mean of the best-fit 
determinations from o ur five transit light cu rves. (3) Derived by fitting the orbital RV curve to out-of-transit RV measurements. 



(4) Derived following Sozzetti et al. (20071. (5) Derived from evolutionary models. (6) Derived by using the expression of the 

mass function. (7) Derived from parameters above. (8) Derived by fitting our RM effect model to in- transit RV measureme nts. 

(9) Derived from spectroscopic analysis. (10) Derived from gyrochronology. (11) Derived following |Cowan fc Agol (20111. 
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